Fast Frictionless Expansion of an Optical Lattice 
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We investigate fast frictionless expansion of an optical lattice with dynamically variable spacing 
(accordion lattice). We design an expansion trajectory that yields a final state identical to the 
initial state up to an irrelevant phase factor. We discuss the effect of additional force and nonlinear 
interaction on the fast frictionless expansion. 



I. INTRODUCTION 



In an adiabatic process, an eigenstate of a time depen- 
dent Hamiltonian follows its instantaneous eigenstate. 
However, adiabatic processes typically need timescales 
longer than lifetimes or coherence times of the system, 
therefore hindering their applicability. This has moti- 
vated researchers to find a way to speed up the process 
to reach the same instantaneous adiabatic eigenstate 
Recently, a new technique has been introduced [2H26l| 
and experimentally realized with ultracold 87 Rb atoms 
in a harmonic trap [ill [l2j . The technique is based on 
engineering the time dependent parameters of the Hamil- 
tonian. The fast schemes to get the same final state as 
the one after adiabatic process_is_ called shortcut to adi 
abaticity @ or superadiabatic 
schemes are not adiabatic in 
In the experiment 



25] 



It is known that fast 
between initial and final 



times. In the experiment [12J, the trap is decompressed 
by a factor of 15 within 35 ms, which is much less than 
typical adiabatic process of a few seconds. Choi et al. dis- 
cussed that optimal cooling of an atomic species may be 
obtained by means of sympathetic cooling with another 
species whose trapping frequency is dynamically changed 
to maintain constancy of the Lewis-Riesenfeld adiabatic 
invariant [l9[ . Recently, nearly perfect fidelity for a two- 
level quantum system was achieved experimentally (25[. 
The fast transitionless mechanism has also been extended 
to the Gross-Pitaevskii (GP) equation [26(. Most of the 
investigations so far have been based on ultracold atoms 
in a harmonic trap. 

In this paper, we will explore the possibility of fast fric- 
tionless expansion for an optical lattice. We will show 
that a correct choice of expansion trajectory allows us to 
get a final adiabatic state in a non-adiabatic way. Our 
method of obtaining fast adiabatic transformation is to 
use a lattice with dynamically variable spacing (accor- 
dion lattice) [27l - l30| . A continuous change of the lat- 
tice periodicity from 0.96 /im to 11.2 fxm in one sec- 
ond was achieved in an early experiment with an optical 
lattice 28]. Recent experiment with ultracold atoms of 
87 Rb trapped in a two-dimensional optical lattice demon- 
strated the variation of the spacing of the lattice in two 
directions from 2.2/xm to 5.5/im in a few milliseconds 
[30T ] using dual axis acousto-optic deflectors. Here, we 
will give the conditions for fast frictionless expansion for 
an accordion lattice. Such accordion lattices are useful 
since the final lattice spacing can be made large enough 



to be resolved in an experiment. Note that, in principle, 
the lattice spacing can be extended to any size smaller 
than the laser beam waist. This allows imaging of the 
atoms at individual sites. 



II. ACCORDION LATTICE 

In contrast to electrons in a crystal lattice, optical lat- 
tices allow us to change parameters, such as potential 
depth and lattice spacing. The potential depth of the 
optical lattice can be modulated by changing the power 
of the laser. The lattice spacing can be changed by chang- 
ing the wavelength of the laser or by changing the relative 
angle between the two laser beams [27h30] ]. Let us con- 
sider ultracold atoms in an optical lattice with variable 
spacing and potential depth. We assume that the trans- 
verse motion of the atoms is frozen (i.e., we are dealing 
with a 1-D problem). The description of a condensate in 
1-D is based on the Hamiltonian 



H =k +mcos [ 2kL W) 



muj 2 {t) 2 



(1) 



where m is the atomic mass, V(t) is the lattice depth, 
is the optical lattice wave number, cu(t) is the time depen- 
dent angular frequency and A(t) is the scale parameter 
describing the expansion of the lattice. The initial and 
final lattice scalings are given by 



A(f = 0) = 1 , A(f = t f ) = c , 



(2) 



where the constant c is the ratio of the final spacing to 
the initial spacing. 

Initially only the periodic potential is present in the sys- 
tem. When the system starts to expand by changing the 
lattice spacing, the parabolic potential is also added. The 
combined presence of the periodic and parabolic poten- 
tials enables us to get fast frictionless expansion. Finally, 
the parabolic potential is turned off at t = tf such that 



u(t = 0) = w(t = t f ) = 



(3) 



The system is initially in the stationary state. To find 
the time evolution of the state for the Hamiltonian ((T|), 
we introduce a transformation on the wave function as 



V(x,t) = exp (~igx 2 /2 + ^— J gdt^j $ 



(x,t) , (4) 
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with a subsequent transformation on the coordinate 
x = zA(t) . 



(5) 



where g(t) is a time dependent function to be determined 
later. Indeed, the scale transformation on the coordinate 
accounts for the expansion of the condensate. Now let us 
substitute these transformations into the corresponding 
Schrodinger equation for the Hamiltonian ([1]) . Note that 
due to the scaling of coordinate, the time derivative op- 
erator transforms as d/dt — >• d/dt — (A/A)z d/dz, where 
dot denotes time derivation. Choosing 



9(t) 



m A 

Ta ' 



(6) 



yields 



h 2 d 2 V 



h— - 

1 dt ~ ~ 2mA 2 dz 2 
where fl 2 is given by 



Vcos(2k L z)$ 



d 2 A 



$ (7) 



(8) 



To sum up, the effects of lattice spacing are equivalent 
to the effective time dependent mass and the effective 
parabolic potential in the stationary frame. The first 
term in the right hand side of © is due to the external 
one while the second one is induced due to the variation 
of the lattice spacing. Indeed, the last formula explains 
the difficulties of keeping atoms trapped in an experiment 
with an accordion lattice [29|, |30| ■ In the absence of ex- 
ternal harmonic trap, lo 2 — 0, the system is left with the 
induced parabolic potential that depends on the acceler- 
d 2 A 

ation, . Since the expansion starts at t — and stops 
dt 2 

at t = t f , the acceleration takes both positive and nega- 
tive values in the interval < t < tf . Hence the repulsive 
potential with an imaginary frequency, O 2 < 0, has the 
effect of pushing the atoms out of the condensate. To 
tackle the problem, the lattice depth was ramped from 
h x 1.5 Hz to h x 37 Hz in 5 ms in the experiment [30| . 
Here we instead suggest decreasing the lattice depth. The 
first reason for this is that the fast frictionless expansion 
can be achieved by engineering the potential depth as 



V{t) 



Vq 
A 2 (t) 



(9) 



Hence the potential depth decreases as the system ex- 
pands. Secondly, we demand that no effective harmonic 
potential acts on the system, VL 2 — 0. In other words, we 
choose uj 2 (t) in such a way that it cancels the effective 
parabolic potential so that there is only periodic poten- 
tial in the system from t — to t = tf. Hence, we choose 



-A 



_ ^ 2 a 
dt 2 



(10) 



This choice of external frequency guarantees that ultra- 
cold atoms stay inside the optical lattice. 
The conditions (|9I10|) reduce the equation to the fol- 
lowing one 



,9$ 
ih— = 

OT 



n 2 d 2 $ 

2m dz 2 



V cos(2h L z)$ 



(11) 



In the last step, we made another transformation on time 
ft 



= f A- 2 (t') dt'. 
Jo 



It is interesting to observe that the Schrodinger equa- 
tion for the Hamiltonian (JJJ can be transformed into the 
another Schrodinger equation (fTTj) with the constant po- 
tential depth. We can view (fTTj) as describing another 
condensate with a different time parameter. 

p 2 

We find a new Hamiltonian H = h Vq cos(2fciz) for 

2m 

the system described by the Hamiltonian (UJ. The final 
state of H' at tt can be made identical up to a global 
phase factor to the final state of the adiabatic evolution 
with H if we impose some conditions on the scale param- 
eter A(t). By virtue of (|4l6p . we conclude that they are 
given by 



A(0)=A(t,)=0, A(0) = A(t / ) = 



(12) 



where dot is the time derivation with respect to the pa- 
rameter t. The acceleration, A, is set to zero at initial and 
final times since the initial and final external frequencies 
are assumed to be zero (|3I10[) . These conditions together 
with the conditions © guarantee the fast transitionless 
evolution. 

Here we will mainly be interested in three different solu- 
tions for A(t). The first one satisfying the relations (|2I12[) 
is the polynomial solution 



A!(i) = ( C -l) 



1 t 3 t 4 t 5 \ 
^T + 10 ^- 15 ^ + 6 d ' (13) 



and the second one is given by 



A 2 (t) = (c-1) 



1 



t 



c-1 t. 



1 . 2mtt 
- — sm( 

2nir tf 



(14) 



where n is a positive integer. For these trajectories, the 
external frequency can be calculated using (fTOf . In FigTJJ 
and Fig-[2l we plot A and corresponding uj 2 . The dashed 
(solid) curve is for Ai (A 2 with n = 3). We assume the 
lattice is expanded by a factor of 2.5. As expected, the 
sign of u 2 (t) changes in time. There are some differences 
in implementing fast frictionless processes in harmonic 
traps and optical lattices. For a harmonically trapped 
system, we need the harmonic potential to achieve the 
expansion of the system. However, expansion for an op- 
tical lattice can be achieved either by changing the wave- 
length of the laser or the relative angle between the two 
laser beams. The external parabolic potential is neces- 
sary for an optical lattice to cancel the induced parabolic 
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A(t) 




0.2 0.4 0.6 0.8 1 t f 

FIG. 1: A(t) versus t/tf for c = 2.5. The dashed (solid) curve 
is for Ai(t) (A a (t) with ra = 3). 




FIG. 2: o; 2 (t) versus i/t/ for c = 2.5. The dashed (solid) curve 
is for Ai(t) (Aa(t) with rc = 3). 

potential, which pushes the atoms out of the condensate 
since it becomes expulsive for certain time intervals. 

One interesting trajectory was already realized in the 
experiment [30| . where the lattice was expanded by a 
factor of 2.5 in a few milliseconds 

A 3 (t) = 1 + 0.75 ^erf{7.14(^- - 0.5)} + 1^ . (15) 

where erf is the Gauss error function. A3 (t) satisfies ([12"]) 
to a very good approximation. Fig-|3]and Fig44]plot A^t) 
and corresponding lo 2 (i) . As can be seen, w 2 is negative 
between < t < i//2, while it is positive otherwise. 
We propose to repeat the experiment [28] with the same 
trajectory but a different potential depth and a different 
frequency. If they change according to (|9I10[) . then fast 
frictionless expansion can be realized. 

So far, we have studied 1-D optical lattice with vari- 
able spacing. The generalization to higher dimensions 
is straightforward. In this section, we will discuss the 
effects of additional force, harmonic potential and scat- 
tering length on the fast frictionless expansion. Let us 
first consider the presence of an additional force, F(t), 
in the Hamiltonian (fl}. The force can be a constant 
gravitational force or a periodically changing one. The 
latter is experimentally realized by generating periodical 
frequency shift between the two laser beams [31j . The 
dynamical responses of the particle to the constant force 
and periodically changing force are the well-known phe- 
nomena of Bloch oscillations [32-34] and dynamical lo- 
calization [33, respectively. The additional force, F(t) 
transforms as F(t)A 3 (t) in the equation (|lip. Hence, 



A(t) 




0.2 0.4 0.6 0.8 1 t f 

FIG. 3: As(t) versus t/tf. The lattice is expanded by a 
factor of 2.5 from 2.2/im to 5.5/im in a few milliseconds. 
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FIG. 4: oj 2 {t) versus t/tf for As(t) 



the picture for Bloch oscillation and dynamical localiza- 
tion change a lot for the trajectories (I13I14I15I) . The fast 
frictionless expansion without destroying the Bloch os- 
cillation and dynamical localization would be possible if 
the force transforms as F/A 3 . 

Suppose now that the system is initially prepared in the 
combined presence of the periodic plus parabolic poten- 
tials. Hence, the initial and final frequencies are given 
by w(0) = to(tf) = ujq. In this case, the condition (fT2"]) is 
still valid but (flU)) is not. The fast frictionless expansion 
for this case is possible if f2 2 is given by Q 2 = w 2 / A 2 in 
the equation ([7|). The specific external frequency, u)(t), 
satisfying this condition can be found easily using the 
equation (|8|). Let us now discuss dipole oscillation. In 
the experiment for 3-D BEC (36|, the trap in the pres- 
ence of the optical lattice was suddenly shifted. It was 
observed that the subsequent oscillations of the BEC are 
undamped if the initial displacement is small, but become 
dissipative if the initial displacement exceeds a critical 
value. Fertig et al. [38[ induced dipole oscillations of 
1-D Bose gases for different values of the axial lattice 
depth. They observed a transition from underdamping 
to overdamping with increasing lattice depth. Due to 
strong quantum fluctuations, mean-field theories fail for 
the description of the system. Consider an experiment 
like the ones in (3(| [37] to study dipole oscillation for 
an accordion lattice. The amplitude of the small dipole 
oscillation increases with A without leading instability 
since the lattice is expanding (|5|). The fast frictionless 
expansion preserves dipole oscillation. 
Let us now discuss the extension of the fast frictionless 
expansion to GP equation. If the scattering length is 
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changed according to a s (t) = a s (0) / A 2 via Feshbach res- 
onance, then adiabatic final state can be obtained non- 
adiabatically. As a result, we say that the fast frictionless 
expansion is not limited to noninteracting case. 
In conclusion, we have shown that the non-adiabatic ex- 
pansion of an optical lattice in few milliseconds may lead 
to adiabatic final states provided that the expansion tra- 
jectory satisfies (|12[) and the external parabolic poten- 
tial is applied to the system. This prediction can be 



tested through experiments similar to the one described 
in [30]. Since the non-adiabatic protocol preserves the 
initial state after expansion, it allows straightforward 
imaging of the optical lattice. We have also discussed 
the effect of additional forces and harmonic trapping. Fi- 
nally, we argue that fast frictionless expansion is possible 
even in the presence of nonlinear interaction provided 
that the scattering length is decreased proportionally to 
A 2 via Feshbach resonance. 
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